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It is argued that experiments on rotating superconductors provide evidence for the existence of 
macroscopic spin currents in superconductors in the absence of applied external fields. Furthermore 
it is shown that the model of hole superconductivity predicts the existence of such currents in all 
superconductors. In addition it is pointed out that spin currents are required within a related 
macroscopic (London-like) electrodynamic description of superconductors recently proposed. The 
spin current arises through an intrinsic spin Hall effect when negative charge is expelled from the 
interior of the metal upon the transition to the superconducting state. 

PACS numbers: 



There is currently great interest in the study of the 
role of the electron spin degree of freedom in transport 
in solids 1] . In particular the possibility of generating a 
pure spin current in the absence of charge currents and 
magnetization has been suggested 0, 13, |j|i ^nd it was re- 
cently reported that it may have been achieved by optical 
techniques @. 

Another area of recent interest is the possible exis- 
tence of novel states of matter that involve spontaneous 
currents, induced by electronic correlations between 
electrons. Early work on superconductivity predicted 
the existence of such macroscopic charge currents 0, 
but this was proven impossible by a so-called 'Bloch's 
theorem' Q. More recently, states with local charge cur- 
rents that spontaneously break time reversal invariance 
have been proposed to describe high temperature cuprate 
superconductors ^ Q ■ 

Combining the above two concepts, we proposed sev- 
eral years ago that the low temperature phase of cer- 
tain metals and in particular Chromium may be a novel 
state of matter induced by electronic correlations ('spin- 
split state') where parity is spontaneously broken but 
time reversal invariance is preserved lOj, in which the 
system carries a macroscopic spin current but no charge 
current |llj. Such spin currents give rise to electric fields 
which should be experimentally detectable [I^, [l3- R-6- 
cently Wu and Zhang discussed a variety of interesting 
possible states of metals arising from spontaneous gener- 
ation of spin-orbit couplingjlJl, one of which is equivalent 
to the proposed spin-split state. 

As pointed out in our discussion of spin-split states |ll|. 
such spin currents should be insensitive to degradation 
by non-magnetic disorder (hence we termed them 'spin 
supercurrents') for the same reason as supercurrents in 
superconductors [iJI- in the presence of weak disorder, 
quasiparticles giving rise to the spin current state should 
be interpreted as time-reversed extended eigenstates of 
the disordered system rather than states with definite 
crystal momentum k. The question then naturally arises 
whether spin currents also occur in superconductors. 

Note that a single Cooper pair (fc |, —k [) carries a 
spin current hut not a charge current, and that the form 




FIG. 1: Qualitative picture of Cooper pairs in a spherical 
superconductor giving rise to spin currents. The arrow per- 
pendicular to the orbit indicates the direction of the electron 
magnetic moment, and the arrow parallel to the orbit indi- 
cates the direction of motion. 
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naturally allows for the existence of a spin current if 
Wki^k) 7^ (u_fc,w_fe). It is then natural to ask whether 
such symmetry breaking will lower the free energy of the 
system. In this paper we show that while this is not 
the case for the attractive Hubbard model, and likely 
also not the case for electron-phonon models of super- 
conductivity, it is the case when superconductivity is 
driven by a kinetic-energy- lowe ring interaction as in the 
model of hole superconductivity [I3, which has been pro- 
posed to describe all superconducting materials [13. In 
the model of hole superconductivity, pairing is driven by 
an off-diagonal matrix element of the Coulomb interac- 
tion; another off-diagonal matrix element of the Coulomb 
interaction was shown to favor the spin-split state[lOj. A 
qualitative picture of the situation envisaged is shown in 
Figure 1. 

Before entering into the microscopic modeling we 
present a macroscopic argument. Consider a rotating 



simply connected superconducting body. A uniform mag- 
netic field ('London field') develops in its interior|13, in 
the absence of any applied external field, given by 
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with TOe and e the free electron mass and charge respec- 
tively, and w the angular velocity of rotation. This has 
been verified experimentally for both conventional p^ 
and high T^ 19] superconductors, and can be 'explained' 
theoreticallV|.2Q| from the London equation 
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and B ~ W X A, with rig the superfluid density. Eq. (4b) 
says that in the interior of the superconductor the super- 
fluid rotates with the same velocity [vg) as the body. In- 
stead, within a London penetration depth of the surface 
the superfluid 'lags' and this relative motion between it 
and the body provides the electric current that generates 
the magnetic field Eq. (2)[l3|. 

However a superfluid electron rotating with angular 
velocity uj requires a centripetal force 
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which is only half of the Lorentz force provided by the 
magnetic field B in Eq. (2)[21|. Hence this situation 
requires a compensating outward electric force on the 
electron to balance the radial force, which would require 
a higher density of negative charge in the interior of the 
superconductor than near the surface. There is no logical 
reason why such a inhomogeneous charge density would 
develop in a rotating superconductor, in particular for 
why the electrons in a rotating normal metal would sud- 
denly move in when the metal is cooled into the super- 
conducting state. 

Instead we propose the following solution of this co- 
nundrum. Assume the speed of electrons of spin a in the 
rotating superfluid is given by 



Va = o-vq + LU X r 



(6) 



for a given spin quantization axis, with vq >> ur. We 
assume furthermore that the vector wq is parallel or an- 
tiparallel to a; x f. The dif ference in the centripetal 
force for a rotating and non-rotating superconductor is 



i^c(w) - Fc(0) = -2meVaUjf ^ -VaX B 
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i.e. precisely the Lorentz force provided by the generated 
London field in the rotating superconductor, Eq. (2). 
Eq. (6) implies that the charge current in the rotating 
superconductor is given by Eq. (4) and in particular is 
zero in the absence of body rotation, and that a spin 
current exists 
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which is independent of the body rotation and in partic- 
ular is non-zero when the body is not rotating. The spin 
current field Eq. (8) is a function of position. 

The existence of such a 'spontaneous' spin current in 
the absence of body rotation and magnetic field requires 
the existence of an electrostatic field E in the interior of 
superconductors, related to vq by Newton's equation 
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Note that the right-hand-side of Eq. (9) is independent of 
(J because the velocity appears squared. We have recently 
proposed a new macroscopic electrodynamic description 
of superconductors |23| which allows for the existence of 
such an electrostatic field, satisfying the equation 
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with Al the London penetration depth, and E^ the elec- 
trostatic field originatig in a uniform positive charge den- 
sity po throughout the volume of the superconductor. 

Conversely we may argue that if an electrostatic field 
exists in the interior of superconductors as predicted by 
the theory of hole superconductivity 23], this requires 
the existence of a velocity field for the supercarriers, and 
consequently requires that each superfluid carrier have 
a 'time-reversed' partner of opposite spin moving with 
opposite velocity so that no charge current results in the 
absence of magnetic fields, and hence leads naturally to 
the Cooper pair concept as well as to the conclusion that 
macroscopic spin currents should exist in the supercon- 
ducting state. 

Next we turn to the microscopic theory. The pairing 
interaction in the model of hole superconductivity in a 
hole representation is given byll5l| 
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where the 'kinetic interaction' term a = At/t originates 
in the correlated hopping term in the Hamiltonian with 
amplitude At[l3, and where t is the single hole hopping 
amplitude and efe is the hole kinetic energy. 

Wc have argued elsewhere that the model of hole su- 
perconductivity predicts that negative charge is expelled 
from the interior of the superconductor towards the sur- 
face in the superconducting state|23. This gives rise to a 
macroscopic outward pointing electric field in the interior 
of the superconductor, and in the presence of such a field 




FIG. 2: Critical temperature versus carrier (hole) concen- 
tration for the model of hole superconductivity (bell-shaped 
curves) and for the attractive Hubbard model (monotonically 
increasing curves). The solid, dashed and dot-dashed lines 
correspond to spin-splitting parameter values 6 = 0, 6 = 0.2 
and b = 0.4 respectively. Parameters in the model of hole 
superconductivity are U = 5, K = 1.65, t = th + nK/8 with 
th = 0.03. In the attractive Hubbard model U — —.15, K=0, 
t = 0.04. Energies are in eV and temperature in K. 



the electronic potential no longer has local inversion sym- 
metry even if the underlying lattice structure does, hence 
the spin-orbit interaction necessarily leads to splitting of 
(k,-k) degeneracy. This then leads necessarily to parity 
breaking and {uk,Vk) 7^ (u„fe,f_fe) in the BCS state. 

For simplicity we consider a two-dimensional or quasi- 
two-dimensional system in a sample of cylindrical sym- 
metry and quantize the spin in the perpendicular direc- 
tion. To allow for the possibility of spin splitting we 
assume at the outset for the kinetic energy of a hole of 
spin a the form 
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where we assume unit lattice spacing. The kinetic energy 
Eq. (12) is time-reversal invariant but breaks parity for 
6^0. In addition to being induced by charge expulsion, 
such a 'spin-split' state is favored by a nearest neighbor 
off-diagonal matrix element of the Coulomb interaction 
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with J > 0, which exists in all solids as discussed in [Tol |. 
Of course it can also arise in a crystal without inversion 
symmetry through ordinary spin-orbit coupling. 

The kinetic energy in Eq. (11) is then Cki or C-ki 
which are identical, and solution to the BCS equations 
proceeds identically to the case in the absence of spin 
splitting, except that in the energy versus k relation the 
form Eg . (12) is to be taken. The gap function is of the 
form [13 
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FIG. 3: Condensation energy dependence on spin splitting 
parameter b in the ground state for the model of hole super- 
conductivity (full lines) and the attractive Hubbard model 
(dashed lines) for two values of the hole density. Parameters 
are the same as in Fig. 2. 



with Am, c constants. Hence Afc 7^ A_fe in general, 
unless Am -^ , A^c ^ which only occurs for a ^ 0, 
U < ior the interaction Eq. (11). 

We show in Figure 2 results for the critical tempera- 
ture versus hole concentration for a case with parameters 
appropriate to the physics of the model of hole supercon- 
ductivity ([/ > 0, a > 0) compared with an attractive 
Hubbard model {U < 0, a — 0). Notably it is seen in Fig. 
2 that spin-splitting enhances the tendency to supercon- 
ductivity when the attraction leading to Cooper pairing 
is of kinetic origin, and suppresses it when the attrac- 
tion is of potential origin as in the attractive Hubbard 
model. It is likely that the same qualitative behavior as in 
the attractive Hubbard model will result for models with 
more extended attractive density-density interactions (eg 
Hubbard model with attractive nearest neighbor interac- 
tions) as well as for conventional electron-phonon interac- 
tion models (Holstein, SSH, Frohlich) where the effective 
electron-electron interaction albeit retarded also leads to 
potential energy lowering. 

In the superconducting state, we find that spin split- 
ting raises the condensation energy in the model of hole 
superconductivity and lowers it in the attractive Hub- 
bard model. Results are shown in Fig. 3 for two values 
of the hole concentration. We find that for the entire 
concentration range where Tc is not zero spin splitting 
increases the condensation energy in the model of hole 
superconductivity and decreases it in the attractive Hub- 
bard model. 

Furthermore in the presence of spin splitting the re- 
lation between chemical potential and carrier concentra- 
tion changes. Figure 4 shows the hole concentration in 
the model of hole superconductivity as a function of spin 
splitting for fixed value of the chemical potential, show- 
ing that spin splitting causes an increase in the hole 
concentration for fixed chemical potential. We interpret 
this result to mean that spin splitting enhances the ten- 
dency for the superconductor to expel electrons from its 
interior|23. This tendency exists for all hole concentra- 
tions. 




FIG. 4: Dependence of carrier concentration on spin splitting 
parameter for the model of hole superconductivity for fixed 
value of the chemical potential. Parameters are the same as 
in Fig. 2. 



Of course the model considered here does not take into 
account the large charging energy that would result if 
the hole concentration changed by the amounts shown 
in Fig. 4. In fact, as electrons move out of the interior 
towards the surface the resulting cost in Coulomb energy 
will stabilize the excess electron concentration near the 
surface at a value corresponding to just a few electrons 
per million atoms, as discussed in ref. 23j . 

Furthermore as charge expulsion takes place the result- 
ing internal electric field that builds up further enhances 
the tendency to spin splitting through the spin-orbit in- 
teraction that occurs between the superfluid electrons in 
macroscopic orbits and the macroscopic electric field that 
builds up, as would occur in a 'giant atom'|2J|. 

Finally, our physical picture provides a simple way to 
understand how the macroscopic spin current builds up 
when the system makes the transition to the supercon- 
ducting state. As negative charge is expelled from the 
interior of the superconductor in the radial direction, the 
interaction of the moving magnetic moment of the elec- 
tron with the positive background gives rise to a force 
in direction perpendicular to the motion whose sign de- 
pends on the spin orientation. Such an effect was dis- 
cussed by us in connection with a disipationless intrin- 
sic mechanism for the spin Hall effect and the anoma- 
lous Hall effect in ferromagnetic metals 25j , and may also 
be predicted by a universal intrinsic spin Hall effect re- 
cently discussed pa. l27[. This force deflects the electron 



in opposite azimuthal directions for up and down spin 
electrons I23 such as to give rise to the macroscopic spin 
current. This is shown schematically in Figure 5. 

In summary we have provided macroscopic and micro- 
scopic arguments in favor of the suggestion that macro- 
scopic spin supercurrents exist in the superconducting 
state of all metals. As already pointed out by BohmQ, 
'Bloch's theorem' does not invalidate such possibilities. 
The direction of flow of these currents is determined by 
the sample geometry through spin-orbit coupling to the 
electrostatic field that results from solution of the macro- 
scopic electrodynamic equations J23. Spin currents fur- 
ther stabilize the superconducting state in the model of 




FIG. 5: Electrons are expelled towards the surface when the 
system goes superconducting, and are deflected by the in- 
teraction of the moving magnetic moment (equivalent to an 
electric dipole moment) with the positive background, giving 
rise to a spin current. The figure also shows the resulting 
equivalent dipole configuration near the surface that arises 
from the macroscopic spin current. 



hole superconductivity and explain the value of the 'Lon- 
don field' in rotating superconductors. When an external 
magnetic field is applied the pre-existent spin currents 
will be modified so that they no longer exactly cancel 
giving rise to a Meissner current that screens the ap- 
plied magnetic field. Experimental detection of macro- 
scopic spin currents in superconductors may be possible 
through spin-resolved neutron scattering experiments as 
discussed in ref. [l^] or angle resolved photoemission ex- 
periments with circularly polarized photons as proposed 
in ref. J3 or optical Q or other experiments|24j| . It has 
also been suggested that the physics discussed here could 
have relevance to aromatic molecules l28l 1291 . 
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